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The DO Collaboration has reported an anomalous charge asymmetry in the 
production of same-sign muon pairs at the Fermilab Tevatron. The magni- 
tude of this effect depends on the subtraction of several backgrounds, the 
most notable of which is due to kaons being misidentified as muons either 
through decays in flight or punch-through. The present authors suggested 
a check on such backgrounds consisting of a tight restriction on the muon 
impact parameter 6, to conflrm that this excess was indeed due to Bf^g) meson 
decays. The DO Collaboration has performed a related check applying trans- 
verse impact parameter (IP) restrictions, whose implications are discussed. 
We study background asymmetry predictions for events involving two muons 
with IP bounds which are complementary to each other. These predictions 
may be used in future measurements of the net charge asymmetry from Bf^g-^ 
decays. 

PACS codes: 12.15.Hh, 12.15.Ji, 13.25.Hw, 14.40.Nd 



I Introduction 

Last year the DO Collaboration reported a charge asymmetry of about 40 times the 
standard model value in the production of same-sign muon pairs at the Fermilab Teva- 
tron [HE]. The magnitude of the effect depended on the subtraction of several back- 
grounds, notably muons due to misidentifled or decaying kaons. The cross section for 
K~ on matter is greater than that for on matter, so that decays in flight pose a 
greater source of muons than K~ decays in flight. We suggested that a tight restriction 
on the muon impact parameter h could test whether this asymmetry was indeed due to 
B meson decays [3J. In an updated version of their analysis, employing a larger data 
sample, the DO Collaboration has implemented this suggestion, conflrming their claim 
for a larger-than-predicted charge asymmetry and providing partial separation of B^ 
and Bg contributions to the charge asymmetry The present paper is devoted to a 
discussion of their results and their relation to our suggestion. 
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In Section [TTl we use numbers of dimuon events corresponding to different choices of 
transverse impact parameter (IP) [1] to estimate the effective value of a related parameter 
b > IP. In Section Hm we show that for maximum values of IP < (50, 80, 120) /xm, scaling 
the kaon background by a factor of ~2/3 leads to a vanishing signal for the dimuon charge 
asymmetry, while the net inclusive muon charge asymmetry is shifted away from zero. 
For minimum values of IP > (50, 80, 120) fim, with much smaller kaon background, 
rescaling of that background does not alter the dimuon and inclusive muon asymmetries 
significantly. Section [IV] discusses predictions of asymmetries for two muons involving a 
minimum IP for one and a maximum IP for the other. We conclude in Section |Vl 

II IMPACT PARAMETERS IP AND b 

In their latest study of the like-sign dimuon charge asymmetry at the Tevatron [1], 
the DO Collaboration investigated the dependence of this asymmetry on the transverse 
impact parameter IP, a quantity equivalent to what we called b± [3]. We advocated 
excluding events with an impact parameter b greater than a chosen value 6o- If the 
average impact parameter in b quark decays is (6), the remaining fraction of dimuons 
from pairs of hadrons containing b quarks will be [1 — exp(— 60/ (b))]"^. Based on a study 
by the CDF Collaboration we estimated (6) = 350 fim. For bo = 100 fim, a value we 
advocated in Ref. |3j, this quantity is then 0.062. Such a sample, thus, should be highly 
depleted of dimuons from b decays. 

Requiring events to involve dimuons with a certain maximum IP (= b±) is less 
stringent than requiring them both to have the same maximum value of b, since b > IP. 
The two are related by 

6= [Ip2 + (6||sin7/.)2]i/2 , (1) 

where b\\ and ip were defined in Ref. P]. It is difficult to impose stringent requirements 
on b\\ [SIE], which is why the DO Collaboration chose instead to restrict the parameter 
IP. With their most stringent constraint, we find their sample of dimuons from b decays 
to be reduced by a factor of 6. 

The DO Collaboration studies samples with both muons having IP either less than or 
greater than the values 50, 80, and 120 /im. We have been provided |6j with the effect 
these cuts have on sample sizes. From this information we are able to extract effective 
values of 60 which indeed exceed IP, by factors of 2.7 to 3.6. In what follows we present 
details of the calculation. 

The numbers of dimuon events both of whose muons have IP above 50, 80, and 120 
fim are shown in Table HI while the numbers of events with IP less than 50, 80, and 120 
lj,m are shown in Table [TTl The number of dimuon events due to b decays is 

N^,^ = N>>F>>R, N<<, = N<<F<<R. (2) 

where the quantities Fss are given in Tables XXI and XXII of Ref. fi], while R = 94% [0] 
is the fraction of the "5*5"' sample coming from b hadron decays. Here we use the 
superscripts >> and << to indicate that both muons have IP either above or below 
the given value. The subscript "S'S"' [HISlll] refers to muons both arising from decays 
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Table I: Number of dimuon events with both muons having IP above 50, 80, and 120 



(/im) (106) (%) (10^) 

50 1.680 85.63±3.74 1.352±0.059 

80 1.152 89.88±5.10 0.973±0.055 

120 0.714 91.79±7.65 0.616±0.051 



Table II: Number of dimuon events with both muons having IP below 50, 80, and 120 
fim. 



jp ]\[« p« ]\[« 

(/im) (10^) (%) (10^) 

50 1.527 43.42±3.75 0.623±0.054 

80 2.174 48.76±2.84 0.996±0.058 

120 2.857 53.66±2.68 1.441±0.072 



of particles at short distances from the interaction point. Muons produced by particles 
traveling long distances before decaying in the detector are labeled "L" . 

The total sample of dimuons due to h hadron decays is N^^^i, = N^^^^ + A^^-^^ + N^^i^, 
where the last term is the contribution from events in which one muon has IP greater 
than the indicated value and the other has IP less than the indicated value. We can 
calculate this term from the first two: 

A^.ir. = 2[iV<<,iV>>,]^/^ (3) 

We should get the same value of N^^^h for each value of IP; the values we obtain from 
the data provided to us are (3.81±0.17, 3.94±0.16, 3.94±0.18) xlO^ events for IP = 
(50,80,120) /im. These are within statistical errors of one another. 

The effective values of Bq may now be calculated from the fractions of dimuon events 
due to b hadron decays with IP greater than a given amount: 

[exp{-bo/{b))r = N>>JN,,,, , (4) 

yielding bo/{b) = (0.518,0.699,0.928) for IP = (50,80,120) /xm. With (6) = 350 /xm 
as estimated in Ref. |3j, this yields bo = (181,245,325) /im, values which exceed the 
corresponding ones of IP by factors of (3.6,3.1,2.7). 

For the most stringent constraint, taking events with IP < 50 /im for both muons, 
the DO collaboration is left with N^^^^/N^^^i, = 0.164 ± 0.015, or a sample of about 
1/6 the size of that employed without imposing bounds on IP values. We note that 
while the asymmetry derived in the latter case, A^i = —0.787 ± 0.172 ± 0.093 is 3.9a 
away from zero [1], a considerably larger asymmetry obtained for IP < 50 /im, = 
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Table III: Contributions of background sources to charge asymmetry in like-sign dimuon 
sample for nominal analysis of Ref . |1] . 



Source 


Asymmetry 


FkAk X 102 


+0.633 ±0.031 


F^A^ X 102 


-0.002 ±0.023 


FpAp X 102 


-0.016 ±0.019 


(2 - Fbkg)A X 102 


-0.212 ±0.030 


Abkg X 102 


±0.402 ± 0.053 


A X 102 


±0.126 ±0.041 


(A - Akg) X 102 


-0.276 ±0.067 



— 2.779±0.674±0.694, is nonzero at 2.9cr. The reason for the larger measured asymmetry 
for muons with low IP may be that it is driven more by Bg decays than by decays [3]. 

Ill RESCALING THE KAON BACKGROUND 

In this Section we examine whether rescaling the kaon charge asymmetry background 
from the value determined in Ref. [4] can lead to a vanishing charge asymmetry for 
same-sign muon pairs originating from decays. We will also discuss the implication 
of such rescaling on the inclusive muon charge asymmetry. 

Using the notations of Ref. [1] , the measured like-sign dimuon asymmetry A is related 
to the charge asymmetry from decays A^i through 

A = {FssCb + FsLCb) All + FkAk ± F^A^ ± F^Ap ± (2 - Fbkg) A 

= {FssC, + FsLCb)A\i + A^^^. (5) 

Here FssCb ± FsiCh is an effective dilution factor depending on fractions of dimuon 
events with two short {S) muons and with one short and one long (L) muon, while 
Fx and A^ {x = K,7c,p) are fractions and asymmetries of muons produced by kaons, 
pions and protons, respectively. The term (2 — Fbkg)A [-Fbkg = Fk ± ± Fp, A = 
(— 0.132±0.019)%] represents a contribution from muon track reconstruction asymmetry. 

We begin by quoting in Table IIIII some relevant entries from Table XII of Ref. [1] 
for like-sign dimuon events, which are the main source of statistical weight supporting 
the claim for a charge asymmetry A^i . In Table IIIII the dominant contribution to Abkg 
is from the charge asymmetry in kaon tracks. The final result for the dimuon charge 
asymmetry is given by (A — Abkg) divided by an effective dilution factor, FssCb + FsLCb = 
0.342 ± 0.028, or 

All = (-0.808 ± 0.202 (stat) ± 0.222 (syst))%. (6) 

Thus a reduction of the kaon background by a factor of 0.56 to 0.356 would be required 
to achieve a vanishing value of A — Abkg or of A^i. 
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Table IV: Contributions of background sources to charge asymmetry in like-sign dimuon 
sample for analysis of Ref. |1] with the restriction IP < (50, 80, 120) fim. 



Source 




Asymmetry 






IP < 50 /im 


IP < 80 /im 


IP < 120 /im 


FkAk X 102 


+1.421 ±0.066 


+ 1.203 + 0.053 


+ 1.047 + 0.051 


F^A^ X 102 


+0.016 ±0.051 


-0.001 + 0.002 


-0.003 + 0.039 


FpAp X 102 


-0.028 ±0.041 


-0.029 + 0.035 


-0.027 + 0.032 


(2 - Fbkg)A X 102 


-0.166 ±0.021 


-0.179 + 0.023 


-0.188 + 0.025 


Abkg X 102 


+ 1.243 + 0.096 


+0.994 + 0.082 


+0.829 + 0.077 


A X 102 


+0.715 + 0.083 


+0.683 + 0.069 


+0.555 + 0.060 


(A - Abkg) X 102 


-0.527 + 0.127 


-0.311 + 0.107 


-0.274 + 0.098 



Table V: Contributions of background sources to charge asymmetry in like-sign dimuon 
sample for analysis of Ref. [5] with the restriction IP > (50, 80, 120) /im. 

Source Asymmetry 

IP > 50 /im IP > 80 /im IP > 120 /im 

FkAk X 102 +0.205 + 0.060 +0.104 + 0.074 +0.113 + 0.087 

F^A^ X 102 +0.001 + 0.007 -0.002 + 0.003 -0.001 + 0.002 

FpApX 102 -0.005 + 0.005 -0.004 + 0.002 -0.002 + 0.002 

(2 - Fbkg)A X 102 -0.244 + 0.035 -0.236 + 0.036 -0.237 + 0.037 

Abkg X 102 -0.043 + 0.071 -0.139 + 0.083 -0.127 + 0.093 

Ax 102 -0.302 + 0.079 -0.386 + 0.094 -0.529 + 0.120 

(A-Abkg)xl02 -0.259 + 0.106 -0.247 + 0.125 -0.402 + 0.152 



Ref. [1] considered a set of variations on their nominal analysis entailing the re- 
strictions IP < (50, 80, 120) /xm and IP > (50, 80, 120) /im. We are grateful to the DO 
Collaboration for sharing the respective entries in Table IIIII corresponding to each of 
these criteria [8J. The results shown in Tables HVl and IVl were obtained by summing 
products such as F^Ak over muon pt bins. 

For comparison, we quote in Tables |VT] and IVIII the dominant contributions F^Ak 
and (2 — Fbkg) A using averaged values of A^ and F^ from Tables VII and XXI, XXII in 
Ref. 15 , respectively. We note the insignificant differences between values of background 
asymmetries in Tables IIVI and I VI I and between those in Table |V] and IVIIi Thus, in 
the next section, where we discuss charge asymmetries in dimuon events with other 
IP constraints, we will use averaged values of fractions and asymmetries rather than 
summing their products over muon bins. 

If we wish to rescale F^Ak by a common factor A, we may parametrize 

A - Abkg = Aa- XFkAk , (7) 
with the values oi Aq = A — [F^A^ + F^Ap + (2 — Fbkg)A] and FkAk given in Table 
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Table VI: Dominant contributions of background sources to charge asymmetry in like- 
sign dimuon sample with the restriction IP < (50, 80, 120) fim. We use average values 
of Ak and F^. 



Source 






Asymmetry 








IP < 50 yum 


IP < 80 fim 


IP < 120 /xm 


FkAk X 


102 


+1.37 ±0.05 


±1.16 ±0.04 


±1.05 ±0.04 


(2 — -Fbkg 


;)A X 102 


-0.158 ±0.023 


-0.173 ±0.025 


-0.181 ±0.026 



Table VII: Dominant contributions of background sources to charge asymmetry in like- 
sign dimuon sample with the restriction IP > (50, 80, 120) /im. We use average values 
of Ak and F^. 



Source 






Asymmetry 








IP > 50 yum 


IP > 80 urn 


IP > 120 jjm 


FrAk X 


102 


±0.31 ±0.09 


±0.23 ±0.13 


±0.22 ±0.20 


(2 — -Fbkg 


;)A X 102 


-0.243 ±0.035 


-0.249 ± 0.036 


-0.252 ±0.037 



Ivml 

The effect of rescaling the contribution of FkAk may be seen in Figs. 1-4. For no 
constraint on IP and for maximum values of IP < (50, 80, 120) /xm, a common choice of 
A ^ 2/3 shifts derived values of A — Abkg to within la of zero. With the exception of IP 
> 120 fim, it is always possible to choose a value of A > 0.66 such that A — A^^g is within 
2a of zero. However, for IP > 120 /xm, with A = 1 the value of A-Abkg = -0.402±0.152 
is 2.6(7 from zero, and quite insensitive to rescaling of F^Ak- 

Effects of rescaling the kaon asymmetry background may also be studied in the 
inclusive muon asymmetry a, which is less sensitive than A to a nonzero value of A^i. 



Table VIII: Parameters in Eq. ([7]) for rescaling contribution of kaon background. The 
nominal parameters [3] involve only a very weak restriction on impact parameter (IP). 



Criterion Aq x 102 FkAk x 102 

Nominal ±0.356 ± 0.059 0.633 ± 0.031 

IP < 50 /im ±0.893 ±0.108 1.421 ± 0.066 

IP<80/im ±0.892 ±0.081 1.203 ± 0.053 

IP<120/im ±0.773 ±0.082 1.047 ±0.051 

IP>50/im -0.054 ±0.087 0.205 ± 0.060 

IP>80/im -0.144 ±0.101 0.104 ±0.074 

IP > 120 /im -0.289 ± 0.126 0.113 ± 0.087 
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Figure 1: Dependence of net dimuon asymmetry A—A^^g on scale factor A for the nominal 
analysis of Ref . [1] . Solid, dashed, dot-dashed, and dotted sloping lines denote central, 
±1(7, ±2(7, and ±3 sigma values, respectively. Values of A to left of vertical line give net 
asymmetries differing from zero (horizontal line) by less than 2a. A — is within 2a 
of zero for A < 0.76. 




Figure 2: Dependence of net dimuon asymmetry A — Ay^g on scale factor A for (left) IP 
< 50 /im; (right) IP > 50 /im. Lines as in Fig. [HA — is within 2a of zero for 
A < (0.80,0.66). 
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Figure 3: Dependence of net dimuon asymmetry A — Ay^g on scale factor A for (left) IP 
< 80 /im; (right) IP > 80 fim. Lines as in Fig. [TJ A — is within 2cr of zero for 
A < (0.92,1.00). 



0.50 



0.25 ■ 



^ 0.00 



-0.25 



0.0 



-0.2 



-0.4 - 



-0.6 



-0.50 — 



0.5 



0.6 0.7 O.B 0.! 

Scale factor A 



1.0 



-0.8 



0.5 



0.6 0.7 O.B O.i 

Scale factor A 



1.0 



Figure 4: Dependence of net dimuon asymmetry A — Ay^g on scale factor A for (left) IP 
< 120 /xm; (right) IP > 120 fim. Lines as in Fig. [T] A — Ay,g is within 2a of zero for 
A < 0.92 (left) but for no positive value of A when IP > 120 /xm. 
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The measured values of a — Obkg for all IP constraints are consistent with zero [UIH] . This 
is expected because the contribution of A^i to this difference is proportional to Cb which is 
much smaller than Cb- (See Eq. ([5]) above and Eq. (7) in Ref. |1] with values of Cb and Cb 
given in this reference.) We checked that rescaling the kaon background contribution to 
a by A ~ 2/3 leaves values of a — abkg consistent with zero within la for IP > 80, 120 /xm 
and within less than 3a for IP > 50 fim. However, for IP < 50, 80, 120 fim this scaling 
moves a — abkg significantly {6a) away from zero. Thus we conclude that a solution 
A^i = cannot be obtained consistently for all IP constraints by introducing a common 
rescaling factor for the background kaon asymmetry. 

IV MAXIMUM & MINIMUM IP FOR TWO /i'S 

So far DO has analyzed dimuon events in which the two muons have IP either less 
or greater than the values 50, 80 and 120 /xm. Other events, which have not yet been 
studied in the DO sample, involve one muon with IP less than 50, 80 or 120 fim and a 
second muon having IP greater than the same value. The numbers of events in the three 
classes, N^^,N^^ and N^^^ , are listed in Table HXl for the three IP values. This Table 
also quotes the calculated total number of dimuon events N^°^'^^ for the three IP values. 
The total numbers are in reasonable agreement with each other. Their deviation by 7% 
relative to the number 6.019 x 10® quoted in Ref. [1] for events with no IP constraint 
may be due to a second order effect [6]. 

Table IX: Number of dimuon events with muons having IP either below or above 50, 80, 
and 120 fim. 



IP 


iV« 
^ Vm 






ATtotal 

^ Vm 






(10«) 




(10«) 


50 


1.527 


1.680 


3.203 


6.410 


80 


2.174 


1.152 


3.165 


6.491 


120 


2.857 


0.714 


2.856 


6.427 



Noting that about half of the total like-sign dimuon events have not yet been studied 
experimentally {N^^ /N^^^^ ~ 1/2), we wish to discuss their expected charge asymme- 
tries. For this discussion we need the fraction of background events F^^{x = KjiTjp). 
Using 

N^,^ = 2{N<<N>>y/\ N<> = 2{N<<N>>YI\ = F^N,, , (8) 

one has 

Fx^"- = iF<< F>>Y^\ ix = K,n,p). (9) 

Thus, fraction of background events from kaons, pions and protons in the class <> under 
consideration may be calculated from corresponding fractions given in Tables XXI and 
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XXII of Ref. [4j for the two classes << and >>. These fractions are listed in Table 1X1 
for IP= 50, 80, 120 fim. The three calculated fractions of kaons for events of the class 
<> are equal within la to the fraction measured by DO for events with no restrictions 
on IP |4j , Fft- X 10^ = 13.78 ± 0.38. We note that kaon fractions are about twice 
larger than the fractions F^'^, while F^^ are about twice smaller than these fractions. 

Table X: Fractions of background events of type << and >> from Ref. |1] and fractions 
of events of type <> given by Eq. for IP = 50, 80, 120 /im. 



Fraction 




IP = 50 urn 


IP = 80 /im 


IP = 120 fim 




X 10" 


-2 


28.03 ± 0.95 


23.79 ± 0.74 


21.49 ± 0.62 


F<< 


X 


10- 


-2 


51.72 ± 3.18 


44.26 ± 2.63 


40.47 ± 2.26 


F« 
p 


X 


10- 


-2 


0.77 ± 0.29 


0.66 ± 0.25 


0.59 ± 0.23 


X 


10- 


-2 


6.31 ± 1.73 


4.79 ± 2.59 


4.48 ± 4.05 


F» 

n 


X 


10- 


-2 


9.51 ± 2.36 


6.39 ± 2.95 


4.43 ±3.95 




X 


10- 


-2 


0.11 ± 0.06 


0.03 ± 0.04 


0.03 ± 0.05 


^ K 


X 


10- 


-2 


13.30 ± 1.84 


10.67 ± 2.89 


9.81 ± 4.44 




X 


10- 


-2 


22.18 ± 2.84 


16.82 ± 3.91 


13.39 ± 5.98 


F<> 


X 


10- 


-2 


0.29 ± 0.10 


0.14 ± 0.10 


0.13 ± 0.11 



Results of asymmetries for background sources are given in Table IXIi The calculated 
background asymmetries, for samples involving one muon with IP larger than 50, 80 or 
120 /im and a second muon with IP smaller than the same value, should be compared 
with corresponding future results by DO. They may be used to extract the net charge 
asymmetry for B(^g^ decays A^i from the raw asymmetry A. 



Table XI: Asymmetries calculated for background events of type <> defined in the text 
for IP = 50, 80, 120 /im. 

Asymmetry IP = 50 /xm IP = 80 fim IP = 120 /im 

F^>Ak X 10^ +0.649 ±0.091 +0.521 ± 0.141 +0.479 ± 0.217 

F<>A^x 10^ -0.007 + 0.018 -0.005 + 0.013 -0.004 + 0.011 

F<>ApX 10^ -0.002 + 0.010 -0.001 + 0.005 -0.001 + 0.005 

(2-F4>)AxlO^ -0.217 + 0.032 -0.228 + 0.033 -0.233 + 0.035 

A4> X 10^ +0.423 + 0.099 +0.287 + 0.145 +0.241 +0.220 



V CONCLUSIONS 

We have examined the relation between the impact parameter (IP) selections performed 
in Ref. |1] and those we suggested in Ref. [3] . We find that the minimum transverse IP of 
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(50, 80, 120) yum considered in Ref. [4| is equivalent to a value of the parameter Bq which 
exceeds IP by factors of 2.7 to 3.6. For the most stringent criterion, taking events with IP 

< 50 yum for both muons the sample of dimuons from B(^s) decays is reduced to about 1/6 
the size of the sample involving no bounds on IP. Fractions of backgrounds from kaons 
and pions in the former sample are each about twice as large as in the latter sample. In 
spite of the considerably smaller sample of signal events with IP < 50 fim and the larger 
background, the statistical significance of A — Abkg = (—0.527 ± 0.127)% measured 
for IP < 50 fim is 4.1cr, the same as measured for dimuons with no IP constraint, 
A - Abkg = (-0.276 ± 0.067)%. (See Tables III and IV.) 

We asked whether a rescaling of the parameter FkAk describing kaon background 
asymmetry could lead to annulment of the claimed charge asymmetry for same-sign 
muon pairs from decays. For selections of maximum impact parameters IP < 
(50,80, 120) fim, a rescaling by a factor of A ~ 2/3 led to reduction of the net charge 
asymmetry to within Icr of zero. However, for minimum impact parameters IP > 
(50, 80, 120) /im, with greatly reduced kaon backgrounds, this was not so, and for IP 
> 120 yum, no positive choice of A reduced the net asymmetry to less than 2a from zero. 
In contrast, while introducing a rescaling factor A ~ 2/3 in inclusive muon samples is 
consistent with A^i = for IP > 50, 80, 120 /xm, it leads to a nonzero asymmetry for IP 

< 50,80,120 yum. 

We calculated background asymmetries for dimuon samples in which one muon has 
a maximum IP while the other has a minimum IP, for the three cases IP = 50, 80, 120 
yum. These calculated asymmetries, expected to be confirmed in future studies by DO, 
may be used for measuring A^i in these samples. 

The fraction of dimuons from kaons was seen to decrease by imposing a minimim 
value for the impact parameter IP. For instance, the background asymmetry from kaons 
was reduced from FkAk = (0.633 ±0.031)% with no IP restriction to FkAk = (0.205 ± 
0.060)% for IP > 50 yum. The corresponding number of like-sign dimuon events decreased 
by about a factor 1/4 (see Table IX), and the significance of a nonzero A — Abkg went 
down from 4. la to 2.4(T (see Table V). In future studies of the same-sign dimuon charge 
asymmetry we thus advocate emphasis on reduction of kaon background by choosing 
a minimum value of impact parameter, even if at the cost of statistics. Such studies 
could, in principle, be performed by other collaborations such as CDF at the Fermilab 
Tevatron and LHCb at the CERN Large Hadron Collider. 
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